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We calculate an analogue of the average phase factor of the staggered fermion determinant at 
imaginary chemical potential. Our results from the lattice agree well with the analytical predictions 
in the microscopic regime for both quenched and phase-quenched QCD. We demonstrate that the 
average phase factor in the microscopic domain is dominated by the lowest-lying Dirac eigenvalues. 



I. INTRODUCTION 

The numerical study of QCD at nonzero chemical po- 
tential [i is hobbled by the complex fermion determi- 
nant, which precludes straightforward Monte Carlo cal- 
culations based on a real measure. Novel approaches to 
the problem by means of Taylor expansion UJ1,J3, 0], 
analytic continuation @, 0| , or reweighting [|[ [J are 
consistent in measuring the response to small chemical 
potential around the temperature of chiral restoration, 
but they disagree at higher values of fi. Moreover, there 
will be fundamental technical difficulties at larger \x and 
larger volumes, and at smaller quark masses and lower 
temperatures: In all these regimes the fluctuations in 
the phase of the fermion determinant will become severe 
(if! . Understanding the behavior of the phase of the 
determinant is thus of crucial importance for the inter- 
pretation of lattice data at nonzero chemical potential 
(see also [Hj]). 

A common laboratory for studying features of the 
finite-density theory is entered by making fj, imaginary 
[HI, [HI]. Here the determinant is real and ordinary Monte 
Carlo methods can be used. On the face of it, it appears 
that this would be useless for studying the phase of the 
determinant since the phase is strictly zero for all con- 
figurations. It is perhaps startling, then, to note that 
one can define and compute an analogue of the average 
phase factor at imaginary /i [l5| . This offers an approach 
to studying the strength of the sign problem occurring at 
real chemical potential by means of simulations at imag- 
inary fi. The advantages of this method are that one 
does not have to deal with the sign problem in order to 
measure its strength, and that one does not have to deal 
with eigenvalues that have spread out into the complex 
plane. 

The phase factor of the fermion determinant at real 
chemical potential is given by 



o 2io _ det(£> + hjq + to) 
det(D + A<7o + m)* 
dct(£> + ^70 + to) 



(1) 



det(Z? — fijQ + to) 
Here we have used the fact that when fi is real, the com- 



plex conjugate determinant is obtained by flipping its 
sign. The average phase factor at real chemical potential 
is then 



2i<?( M )\ = / det ( D + A*7o + m) 
I \ dct(Z? — /170 + to) 



(2) 



As in [151 ] we define the average phase factor at imaginary 
/i by simply substituting ifi for [i 



i0(iu)\ = / dct(J> + z/x7o +to) 
\ det(_D — i/i7o + to) 



(3) 



where the parameter /i is still real. Both determinants in 
Eq. ([3]) are now real. We present in this note a numer- 
ical study of the phase factor defined in this manner at 
imaginary /i. In particular we study the dependence on 
the chemical potental, the quark mass, and the volume. 
We work in the microscopic domain of QCD (also called 
the e-regime) 0, [ItJ , where the quark mass to and the 
chemical potential are chosen to fulfill 



\n\ 2 F% < and m£ < 



(4) 



while the Euclidean volume V = L 3 /T (where T is the 
temperature) is taken large, 



VA A Q CD » 1. 



(5) 



Here £ = (ipip) is the chiral condensate and is the 
pion decay constant. 1 Formulas for the average phase 
factor have recently been derived in this regime [la . [l8| . 

The values of S and F K on the lattices we use were 
calculated from a two-point spectral correlation func- 
tion in Refs. and [2(j. Given these values, we can 
make parameter-free comparisons of the numerical mea- 
surements of the average phase factor with the analytical 



1 In fact we take \u\ 2 F% ~ 1/V = T/L 3 , which implies that fj,/T < 
1 when L is large. Thus the periodicity in the imaginary chemical 



potential noted by Roberge and Weiss Il3t will never appear 
the microscopic domain. 



in 



2 



predictions. In all cases we study the agreement is within 
the statistical errors. This confirmation of the analytic 
predictions at imaginary fi gives support to the analytic 
results derived for real fj, where direct lattice tests are 
harder to obtain. Our results also show that the average 
phase factor in the microscopic domain is dominated by 
the lowest-lying eigenvalues of the Dirac operator. 

Below we will recall the theoretical predictions in the 
microscopic domain and their analytic properties as func- 
tions of the chemical potential. Then we will compare 
these to the lattice data. 



fermions given an imaginary isospin chemical potential. 
Here we present the microscopic formulas obtained for 
these two cases. 

Below we work with standard unimproved staggered 
fermions at a coupling where the chiral symmetry break- 
ing pattern is identical to that of the continuum theory 
in the sector of zero topological charge (see e.g. [U). We 
therefore present the analytical predictions in the trivial 
topological sector. 



II. ANALYTICAL FORMULAS 



A. Quenched theory 



Our numerical results are based on data obtained in a 
quenched ensemble and in an ensemble with dynamical 

I 



The average phase factor at real chemical potential in 
the quenched theory is given by [HI, Il8j 



= 1 -4/rI (TO)i£o (m) 

/ JV/=0 



4/7 



r,2 



dxxe sa 5 K yj~2 ) [h{x)mh(m) - xh{x)I {rh)] , (6) 
I 



where we have defined microscopic variables via 
fh = mY,V and fl — ^iF^VV. 



(7) 



The formula ((6|) cannot be continued to imaginary /i be- 
cause of the essential singularity of the last term at /x = 0. 
The non-analyticity has its origin in the inverse determi- 
nant of the non-hermitian Dirac operator in Eq. (|2|) — it 
is due to the ei genv alues with real part larger than the 
quark mass [l5l \22\ . For purely imaginary /i the eigen- 
values remain on the imaginary axis and are thus always 
inside the quark mass. The nonanalytic term therefore is 
not expected to appear at imaginary fj,. Indeed, a direct 
calculation at imaginary /x gives [15| 



N f =0 



= l + 4fi 2 Io(m)K (m). 



(8) 



As expected the analytic continuation of this result from 
imaginary values of the chemical potential back to real 
values gives just the first two terms in Eq. ©. For real 
/i the eigenvalue density outside the quark mass is highly 
suppressed if /i < Consequently the first two 

terms in ^ are dominant when fi < (see Fig. [1}. 

The measurement of the average phase factor at imagi- 
nary values of fj, is therefore indicative of the strength of 
the sign problem for real /x < m T /2. 



Two dynamical flavors with isospin chemical 
potential 
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FIG. 1: The quenched average phase factor as a function of 
(2fj,/m w ) 2 for m = 3. The average phase factor from fj, 2 < 0, 
Eq. ([H, gives the solid line when continued via i\i — > /i. This 
line agrees with the direct result for real /x, Eq. JB), when 
(i < m,r/2. For larger values of p the essential singularity in 
Eq. ([6]) at fi = becomes important and the analytic contin- 
uation fails to reproduce the true result. The same occurs in 
unquenched QCD. 



In this theory, also known as the phase-quenched the- 
ory, there are two flavors of fermion with opposite val- 



ues of chemical potential ±/i. The average phase of the 
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single-flavor determinant is given (for real fx) by 



i 2WM ) 1 lt = (9) 
(|det(C + M7o + m)| 2 ) 



where the averages on the right-hand side are quenched 
averages. This simplifies to 



(det(D + fij + to) 2 ) _ Z Nf=2 



(|det(£) + /i7o + m)| 2 



Z 



1+1* 



a ratio of partition functions in ensembles with dynamical 
fermions. These can be evaluated exactly in the micro- 
scopic regime [13, to obtain 



,2i0(AO\ 



7 2 (m)-/ 2 (m) 



/i+i* 2e 2 A 2 J^dtte- 2 P- H2 I^(rht) 



(10) 



for real fi, with scaling variables as defined in Eq. Q 
(note that the numerator, which is Zjv,=2> does not de- 
pend on fi). In this case the average phase factor is an- 
alytic at fi = 0, because the inverse determinant in the 
numerator of Eq. {9} has been canceled. The average 
phase factor therefore is 



,2ifl(»M)\ 



J 2 (m)-J 2 (m) 



(11) 



/i+i* 2e" 2 A 2 f^dtte^lgimt) 
for imaginary chemical potential «/^. 

III. NUMERICAL DATA 

In the lattice theory, we work with standard, unim- 
proved staggered fermions and introduce the chemical 
potential using the Hasenfratz-Karsch prescription [24| . 
As in the continuum, each operator D ± ifijo is anti- 
hermitian, and anticommutes with 75. Each operator's 
eigenvalues therefore come in pairs of opposite sign on 
the imaginary axis. For each gauge field configuration 
we measure the two sets of eigenvalues defined by 



[D{A)+ifx l0 ]ipf =i\ { i ) i,f, 
[D(A) - ifijo]ip U) = iA^'V-*- 



(12) 
(13) 



where the products run over positive A's only. The equal- 
ity of ($_| ) and ($ |_) follows from charge conjugation 

symmetry, which exchanges A + with A_ and is true in 
both the quenched and the phase-quenched theory (but 
not in the unquenched theory). With finite statistics, the 

two estimators 4>_| and $ |_ give distinct measurements 

of the average phase. 



A. Quenched theory 

Our quenched eigenvalue data are a combination of 
data used in Ref. [19| with new data. We simulated the 
SU(3) gauge theory with the plaquette action at (3 = 5.7 
on lattices with 8 4 and 12 4 sites. For each gauge con- 
figuration we calculated the smallest 24 pairs of positive 
eigenvalues A± . 

On the smaller lattice, we chose two values of the chem- 
ical potential: afi — 0.01 and afi = 0.1. In Ref. [l|| we 
determined F v and S, giving the scaled values fi = 0.159 
and fi = 1.59, respectively; also T,V/a — 1039, which per- 
mits converting the lattice masses am into scaled masses 
to. With these values we can make a parameter-free com- 
parison between the analytical prediction ([5]) and the 
measurement of the (truncated) average phase factor. 
The result for afi = 0.01, based on data from Ref. [Tgj | . 
is shown in Fig. O The agreement is well within the 
statistical error bars. 




0.8 1 1.2 



Thus det(Z) ± ifijo + to) = IXj-(iA± + m )- Combining 
positive with negative A's, the average phase factor is 
given by 




A°' )2 +to 2 ' 



w A^' )2 + to 2 



n 



(*+-) (14) 
<*-+> (15) 



FIG. 2: The average phase factor calculated via Eq. (|14[) 
(circles) and via Eq. (|15p (crosses) from the lowest 24 pos- 
itive eigenvalues of the staggered Dirac operator on an 8 4 
lattice, for fixed chemical potential afi = 0.01 as a function 
of the scaled quark mass; 4000 configurations. The full line 
is the analytical microscopic prediction using the parameters 
determined in Ref. 

The value afi = 0.01 was chosen in Ref. [l!| to en- 
sure that the scaling variable ft, was roughly 1/10, since 



4 



this facilitates the determination of F^. At such small 
values of fi the shift of the eigenvalues due to the imagi- 
nary chemical potential is small compared to the average 
level spacing. 2 This was our motivation in running a 
new quenched simulation in order to calculate eigenval- 
ues with the larger value afi = 0.1. The result is shown 
in Fig.(3J again the agreement with Eq. ([5]) is good. Note 
that in this case /J.L = 0.8 which is pushing the limit of 
the microscopic domain. 



1.03 
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FIG. 3: As in Fig. [2]but with afi = 0.1; 



configurations. 



Since m is not a dynamical parameter of the quenched 
simulation, the data points for different values of m are 
highly correlated, and the variation among them is sys- 
tematic, not statistical. In effect, m determines which 
eigenvalues contribute the most to the average phase in 
Eqs. (|14[) and (|15[) . When m is small, all eigenvalues 
contribute to the phase; the smallest eigenvalues cause 
large fluctuations, while the large eigenvalues show but 
little dependence on fi and thus contribute little to the 
ratio. When m is large it eliminates the effect of small 
eigenvalues on the phase. Therefore as m approaches the 
largest of the calculated eigenvalues there is a system- 
atic deviation from the analytic curve. We illustrate this 
by cutting the number of calculated eigenvalues to 10. 
Figure U shows the deviation from the analytic curve as 
m£V approaches A^EV ~ 30. 

For the larger lattice, with 12 4 sites, we use data from 
Ref. with a/i = 0.002. In Ref. pi] we determined 
that fi = fiF^VV = 0.159 and EV/a = 5259. Again we 
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N=0 8 4 \iFV" 2 =0.1S9 LV=1039 lOev's 
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FIG. 4: Same eigenvalue data as in Fig.[2j but retaining only 
10 eigenvalues. As the quark mass approaches the largest 
retained eigenvalue (A' 10 ^£F ~ 30) we observe that the error 
band drops below the analytic curve (thin line). 



can compare the analytical prediction 1(5)) to the measure- 
ment of the truncated average phase factor without free 
parameters, as shown in Fig. [5j The agreement is clear, 
although statistically weaker than for the smaller lattice. 



o <<E> + > 
X <<£ > 



Nj=0 W \iFV" 2 =0.0717 ZV=5259 




FIG. 5: As in Fig. but for the 12 4 lattice; 2000 configura- 
tions. 



2 For a real /i of this order the width of the strip of eigenvalues in 
the complex plane is much smaller than the average level spacing 



along the imaginary axis 1231 125| . |2 ! 



B. Two dynamical flavors 

For the two-flavor theory with imaginary isospin chem- 
ical potential we use eigenvalue data calculated in the 
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long runs of Ref. [2(| , which should be consulted for nu- 
merical details. The lattice theory is based on the pla- 
quette action and unimproved staggered fermions, sim- 
ulated at f3 — 4.2 with volume 6 4 . We have data for 
two masses, am = 0.002 and 0.005, and the imaginary 
chemical potential coupled to isospin is afi = 0.0055. A 
determination [20J of aF^ gives /t = 0.1338. 

For the smaller quark mass, am — 0.002, we calcu- 
lated |20[] that to = mYy = 3.318. From the lowest 24 
eigenvalues we measure 

= 1.0092 ± 0.0020 (16) 
= 1.0082 ±0.0020. (17) 

The analytical prediction for the average phase factor 
CD} is 

(e 2l6 ^ l) ) = 1.0094, (18) 
\ / 1+1* 

in agreement with the numerical averages. 

For the larger quark mass (ma = 0.005) we have 
to£V = 8.295. From the lowest 24 eigenvalues we com- 
pute 

= 1.0036 ±0.0030 (19) 
(*-+)i+i* = 1-0032 ±0.0030. (20) 

The analytical formula (fTTj) gives 

(e 2ie ^A = 1.0041. (21) 
\ / 1+1* 



Here also we see satisfactory agreement with theory, al- 
though the deviation from unity is not really significant. 

IV. CONCLUSIONS 

We have computed an analogue of the complex phase 
factor of the fermion determinant at imaginary values 
of the chemical potential. The results have been used 
as a parameter-free test of the predictions from the mi- 
croscopic domain of QCD. In the quenched as well as in 
the dynamical cases studied the numerical and analytical 
results are in agreement. Furthermore, we have demon- 
strated that the phase factor in the microscopic domain 
is dominated by the low-lying Dirac eigenvalues. 

We hope that this first numerical study of the phase 
factor at imaginary chemical potential will encourage 
studies also outside the microscopic domain. 
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